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Highlights 


© Transient variation during radon exhalation rate measurements was simulated by 


CFD. 
@ The effects of porous medium and measurement parameters were examined. 


@ The discrepancy between various radon exhalation rates was quantified and 


analyzed. 
@ Simulation results are expected to optimize the closed-loop measurement method. 


Abstract: Accurate measurements of the radon exhalation rate helps identify and 
evaluate radon risk regions in the environment. Among these measurement methods, 
the closed-loop method is frequently used. However, traditional experiments are 
insufficient or cannot analyze the radon migration and exhalation patterns at the gas— 
solid interface in the accumulation chamber. The CFD-based technique was applied to 
predict the radon concentration distribution in a limited space, allowing radon 
accumulation and exhalation inside the chamber intuitively and visually. In this study, 
three radon exhalation rates were defined and two structural ventilation tubes were 
designed for the chamber. The consistency of the simulated results with the variation 
in the radon exhalation rate in a previous experiment or analytical solution was 
verified. The effects of the vent tube structure and flow rate on the radon uniformity in 
the chamber; permeability, insertion depth, and flow rate on the radon exhalation rate; 
and the effective diffusion coefficient on back diffusion were investigated. Based on 
the results, increasing the insertion depth from 1 to 5 cm decreased the effective decay 
constant by 19.55%, whereas the curve-fitted radon exhalation rate decreased (lower 
than the initial value) as the deviation from the initial value increased by 
approximately 7%. Increasing the effective diffusion coefficient from 2.771077 to 
7.77*10°° m? s! made the deviation expand from 2.14% to 15.96%. The conclusion is 
that an increased insertion depth helps reduce leakage in the chamber, subject to 
notable back-diffusion, and that the closed-loop method is reasonably used for porous 
media with a low effective diffusion coefficient in view of the back-diffusion effect. 
The CFD-based simulation is expected to provide guidance for the optimization of the 
radon exhalation rate measurement method and, thus, the accurate measurement of the 


radon exhalation rate. 
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1 Introduction 


Radon, a colorless and odorless naturally radioactive gas, is the second leading 
cause of lung cancer.!!: 7] Widely found in building materials, soil, rocks, and uranium 
mill tailings,!*! concerns have been raised regarding radon released from these media 
into the ambient air. An accurate radon exhalation rate can be used in environmental 
assessments to quantify and identify radon-risk areas.!°! 


Currently, radon exhalation has been studied,'”"'7! and the methods commonly 
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used for radon exhalation rate measurement include activated carbon,| l open 


15, 16] and closed loop." For a closed-loop system, the effects of leakage from 


loop,|! 
the accumulation chamber, back-diffusion at the gas—solid interface, and sampling 
flow rate during radon exhalation rate measurement are nonnegligible.!'® 7° The 
increased insertion depth proved to be effective in reducing chamber leakage," while 
ignoring the back-diffusion effect. The back-diffusion phenomenon is caused by the 
radon returning to the porous medium as the radon level increases in the accumulation 
chamber.!*?] This affects radon exhalation on the surface of the medium, as well as the 
measurement of the radon exhalation rate.*?7! Moreover, convection (or advection) 
driven by the pressure gradient is an important factor influencing radon exhalation.?* 
2] Conventionally, radon exhalation rate measurements for porous media with 
unknown physical parameters do not show a discrepancy between the measured and 
initial values. The limitations of experimental tests, such as the control of 
environmental factors and experimental parameters, make it difficult or time- 
consuming to perform the experiments. Although the analytical solutions of 1D and 
2D models can reflect radon migration patterns in steady or transient states, they are 
subject to various assumptions such as ignoring the insertion depth, lateral diffusion, 


specific dimensions of the accumulation chamber, or porous media, and constant 


uniform concentration distribution inside the chamber. !3%! 


Diffusion driven by the concentration gradient and convection driven by the 
pressure gradient are responsible for radon migration in porous media. However, the 
traditional measurement method cannot provide details of radon migration in the 
chamber and gas flow in the porous medium. These are the key data for analyzing the 
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underlying reasons for the deviation of the measured radon exhalation rate from the 
true value. The advancement of computational fluid dynamics (CFD) has made it 
possible to numerically solve the Navier-Stokes equations and has been maturely 
applied to scenarios such as soil, atmosphere, indoor environments, underground 


34°38] providing researchers with 


space pollutant migration, and nuclear reactor system, |! 
solutions to complex scientific problems. Relevant simulations have been performed 
in previous studies to predict the radon concentration distribution and migration in a 
limited space,?*3*! indicating good consistency between the experimental and 
simulation results. In addition, a two-phase flow and heat transfer simulation with a 
simplified porous media model was established for engineering purposes, and the 
results showed good agreement with the experimental data.'*"! Therefore, the CFD- 
based 3D numerical simulation exhibited spatial and temporal radon migrations in the 


accumulation chamber and porous media and revealed the variation in radon 


exhalation at the gas—solid interface during the measurement. 


In this study, two types of vent tube structures (open-ended and half-open) in the 
accumulation chamber were modeled, and three radon exhalation rates ( initial, 
transient, and fitted radon exhalation rate) were defined to quantify the radon 
exhalation behavior. Using uranium mill tailings as the research object, a geometric 
model of a radon exhalation rate-measuring device and a mathematical model for 
radon migration were established. The CFD-based technique was employed to 
demonstrate the radon migration patterns under the coupling of multiple factors (flow 
rate, insertion depth, diffusion coefficient, and permeability), and thus analyze the 
variations of the fitted and transient radon exhalation rates at the gas—solid interface. 
The simulated results are consistent with previous experimental phenomena and are 
expected to optimize the design of an accurate measurement of the radon exhalation 


rate using the closed-loop method. 


2 Materials and methods 
2.1 Theory of radon diffusion and exhalation in porous media 


One-dimensional steady-state diffusive migration of radon in a porous medium 


can be described using the following equation: 


ge 
p45, Ce, (1) 
Ox 7] 
where C is the radon concentration in the air, Bq m”; x is the vertical depth of the 
porous medium, m; D is the diffusion coefficient of radon in the porous medium, m? s- 
1; 7 is the porosity, %; a is the free radon production rate, Bq m? st; and Arn is the 


decay constant of radon, 2.1<10° s. 
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Fig. 1. One-dimensional model of radon diffusion in porous medium 
The boundary conditions in Eq. are C(0)=C4 and 0C/ ôx) =0 at the top and 
bottom of the porous medium, respectively (Fig. 1). The radon concentration in the 
boundary air was significantly lower than that in the porous medium, making it 
possible to assume C,=0 Bq m°. Upon substituting of JD/A for L, the analytical 


solution of Eq. is: 
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The radon exhalation rate for diffusion is then theoretically deduced: 


w <a, (3) 
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C(x) = re ( 


where Ez is diffusive radon exhalation rate, Bq m? s“. 


2.2 The closed-loop method and device for radon exhalation rate 
measurement 

The closed-loop method is commonly used in the measurement of radon 
exhalation rate, which involves the accumulation of radon in the chamber (Fig. 2). 
Radon exhaled from the surface of the porous medium flowed into a radon detector 
(e.g., RAD7 from Durridge Company Inc., USA) driven by an air pump and returned 
to the chamber through the flow-in vent tube. 
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Fig. 2. Connection diagram of the system. 


Inside the accumulation chamber, the accumulation of radon is described as: 


dC ES 
— =—*-), C-A,C-AC 4 
dt V ‘Rn A, 1 ( ) 
Where ¢ is the accumulation time, s; E is the radon exhalation rate, Bq m? s!; Sa 
is the exhaled area of the porous medium, m7; /, is the back-diffusion rate of radon in 
the chamber, st; A; is the leakage rate of radon in the chamber, sl; V is the volume of 


the chamber, m?.j. is the effective decay constant, a sum of radon decay, back- 


diffusion, and leakage rate (AA=ArntAsti), S1. 
The boundary conditions in Eq. , C(0)=Co. Upon substitution of Co for 0 Bq m? 
in approximation, the curve fitting solution for Eq. is 


ES “at 
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The linear fitting solution for Eq. isl! 


ES, 


C(t)=—4t, 6 
Da (6) 


In subsequent calculations, the fitted radon exhalation rates of Efex and Epin 
were obtained by curve fitting (Eq. ) and linear fitting (Eq. ), respectively, with the 
criterion R* > 0.9: The calculation of R? involves the values of the residual sum of 
squares (RSS) and the total sum of squares (TSS), the equations of which are as 


follows: 
R=1-— =1-____, (7) 


where n is the total sum of the data points, y; is the actual data point, f; is the 


value obtained from the fitted curve, and y is the mean value of all data points. 


3 Methods of numerical simulations 
3.1 Measuring device and geometric model 

The three-dimensional geometry (Fig. 3) for the model of the accumulation 
chamber (AC) had internal dimensions of 30 cmx30 cmx10 cm with a thickness of 4 
mm, inserted into uranium mill tailings (a height of 50 cm filled in a container of 100 
cmx100 cmx60 cm). Two vent tubes (06 mm) were designed on the diagonal of the 
chamber, with 3 cm from each hole to the inner wall of the chamber. Besides the AC 
model, a reference (R) model, a container dimension of 100 cmx100 cmx60 cm, was 
established to determine the reference radon exhalation rate by the closed-loop 
method for uranium mill tailings (a height of 50 cm) filled in the container being 


sealed but with two holes (06 mm, 50 cm apart) on the top cover. 
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Fig. 3. Three-dimensional geometric model of the accumulation chamber from (a) top (b) lateral, and (c) overall 
views. 


Considering the disturbance of the porous medium caused by the jet from the 
inlet vent tube, two vent tube schemes (open-ended and half-open) were designed. 
The first type is a vent tube with an open end. The other was sealed at the bottom end 
of the vent tube with four holes (2mm, 2 cm apart, and a normal direction of 90° of 
two adjacent holes) open on the tube wall to avoid the direct impact of the jet from the 


vent tube on the porous medium. 


The chamber had a thickness of 4 mm, and its volume was subtracted from that 


of the porous medium upon insertion. 


The geometry was discretized into polyhedral and hexahedral elements using 
Fluent Meshing to reduce memory and computing time. The mesh was refined in the 
region (e.g., the gas—solid interface and local opening of the vent tube) where rapid 


changes in the radon concentration and air flow were expected, as shown in Fig. 4. 
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Fig. 4 (a) Discrete mesh from overall view. (b) The sectional plane of discrete mesh where vent tubes are included. 


3.2 Governing equations 
The gas flow, which is regarded as incompressible, in the chamber and porous 
medium followed the laws of conservation of mass and momentum. Following the 
law of conservation of mass, the mass conservation equation is described by!*! 
V-v=0, (8) 
where ¥ is the physical velocity vector, m s^. 


Following the law of momentum conservation, the momentum conservation 


equation is described byl“! 


Be? 4 9.(p,58)=-Vp + p8 +y- E-F, O) 


where pa is air density, kg m’; p is the static pressure, Pa; g is the gravitational 
vector, m s?; and F , the source term of external body forces, N m?, describes the 


viscosity of the porous medium, and is only present in the porous zone. The laminar 


flow is described as follows: 


; (10) 


where K is permeability, m?; u is the dynamic viscosity, Pa s; 7 is the porosity of 


medium. 


The viscous stress tensor (T ) is described by!*5! 


T=(u+m)(V5 +V"), (11) 


2 
where ju is the turbulent viscosity, Pa s, defined by 4, = pC, a (Ca is a constant; 
E 


k is the turbulent kinetic energy, m? s?; and e is the turbulent dissipation rate, m? s°). 


The governing equation for radon migration is expressed as 


oC OC C Sr 

— =D- -v t, (12) 

Ot Ox Ox, n 

v, si yp, (13) 
nu 


where D is the radon diffusion coefficient (m? s!); v; corresponds to the physical 
velocity of convection in the x, y, and z directions (m s`); Srn is the source term (Bq 


m° s'); P; corresponds to the pressure in the x, y, and z directions (Pa). 


In a porous medium, Sz, is determined by the source and decay terms 


(S,, =-AnC,, +@), and D=Dm. In the accumulation chamber, the porosity is set to 1; 
Srn is determined by the decay term (Sz, =—A7C;,, ) and D=Dm+D;(D; is the turbulent 


diffusion coefficient, m? s', defined by where turbulent Schmidt number (Scr) 


psc, 


is 0.7 6l; Dm is molecular diffusion coefficient of radon in air, 1.05x105 m? s1). 


3.3 Parameters for simulation 

The simulated scenarios involved reference (R) and accumulation chamber (AC) 
models. In the AC model, the effects of the vent tube structure and permeability of the 
porous medium on radon exhalation were analyzed. Furthermore, the effects of 
insertion depth (H1/H3/H5 corresponding to insertion depths of 1/3/5 cm), flow rate 
(0.5/1/2 L min’), and effective diffusion coefficient (D1/D2/D3/D4 corresponding to 
2.77*10°7/7.77X10-7/2.77*10°°/7.77x10° m? s) were also tested. The effective 


diffusion coefficient of radon in a porous medium is defined asl“?! 
D,=nD=ntD,,, (14) 


where z is the tortuosity factor, 0.66 (dimensionless). 47] 


In the R model, the radon exhalation rate obtained from the reference container is 
used to evaluate the rate measured using the AC model. The design of the R model 
helps improve the reliability of the simulation, considering that the reference radon 
exhalation rate is commonly used in laboratory tests to evaluate the radon exhalation 
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rate obtained from the AC model. 


The free radon production rate (a) is defined byl“?! 
& = AP, Ára Ern» (15) 


where ps is the density of the porous medium, kg m°; Aza is the activity 
concentration of radium, Bq kg!; and Ern is the radon emanation coefficient, 


dimensionless. 


Radium activity concentration of uranium mill tailings commonly ranges from 
less than 5 kBq kg! to as high as 10 kBq kg"!,/**! with the emanation coefficient 
ranging from 0.1 to 0.35, dry-bulk density of 1800 kg m°, and averaged porosity of 


0.4 measured using the drainage method for samples in the laboratory. 


Laminar flow in a porous medium follows Darcy’s law, and permeability, which 
is an intrinsic property of a porous medium, is defined by the Kozeny-Carman 
empirical equation!*”! 


D 


P 


~ 180 0-7)” 


(16) 


where D, is the average grain size (mm) and K is the permeability, m7. 


After sieving the tailings sand in the laboratory, Dp was determined to be 0.5 mm 
and K was calculated to be 2.47x10°'° m’. The physical parameters for the simulations 


were determined, as listed in Table 1. 


Table 1 Physical parameters used in the simulation 


Physical parameters Value 

Free radon production rate, Bq m! s! 3 112,50] 

Porosity 0.4 BY 

Effective diffusion coefficient for 2.77x107/7.77x10; 

radon, mê s“! 112.77 10°/7.77 106 U9! 
Permeability, m? 1x10°/1 x 1071/1 x 10-1 [52 

Flow rate, L min”! 0.5/1/2 

Insertion depth, cm 1/3/5 


The finite volume method (FVM) was adopted for numerical calculations. The k- 
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e turbulence model was set for the chamber zone due to confined jet from the four 
openings on the surface of vent tube, whereas the flow in the porous media was 
laminar model. The empirical constants in the k-e turbulence model equations!**! are 


as follows: C= 0.09, Cie = 1.44, Coe= 1.92, ox = 1.00, and os = 1.30. 


To simplify the model, the following assumptions were made: the leakage of the 
container was not considered, and the porous medium was homogeneous, with air as 


an incompressible gas, ignoring the impact of gravity and temperature. 


3.4 Initial value and boundary settings 
The radon concentration distribution of the porous medium in the steady state 
was calculated before simulating the radon accumulation in the measurement, in 
which the radon concentration in the air was set to 0 Bq m*. Transient mode was 
adopted upon the measurement, with initial radon concentration at 0 s in the chamber 


set at 0 Bq m®. The detailed boundary settings are listed in Table 2. 


Table 2 Boundary settings in the simulation. 


Boundary Settings 
Inlet/Outlet Velocity- Inlet: 0.5/1/2 L min`, normal to boundary. 
Inlet radon concentration equals to that from 
outlet 
Outlet: -0.5/1/2 L min“, normal to 
boundary. 
radon concentration gradient equals to 0 Bq 
m4 
The surface of porous Pressure- 0 Pa; radon concentration equals to 0 Bq m 
medium not covered by Inlet 3 
the chamber 
Openings of the vent Interior — 
tube inside the chamber 
Gas-solid interface Interior — 
Wall of the geometry Wall Radon flux equals to 0 Bq m? s“! 


3.5 Mesh independence 

The three-dimensional geometric model was discretized in coarse, medium, and 
refine (0.5 M, 1.2 M, and 2.3 M cells, respectively) mesh and the radon exhalation 
rate for each mesh was further calculated as 1.4676593, 1.4674347, 1.4674092 Bq m? 
sl, respectively, shown in Fig. 5. The decline in the radon exhalation rate slowed with 
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an increase in the sum of cells. The criterion for discretizing the mesh is to minimize 
the error between the numerical and analytical results while maintaining the 


calculation at a lower cost. Thus, a medium mesh was used in subsequent simulations. 


The calculation is considered converged with the residual less than 10°!° and the 


differential value of radon exhalation rate of two consecutive iterations less than 107. 
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Fig. 5 Mesh sensibility analysis. 


4 Data processing 


As the FVM method was adopted in the numerical calculation, the volume- 
averaged radon concentration (C\-ag) in the chamber was obtained using discrete data 


in the approximation, as follows: 


(i an (17) 
V, 


i=l 


where n is the total number of cells inside the chamber, V; is the volume with a 


cell index i, m*, C; corresponds to the radon concentration in cell Vi. 
The area-averaged radon concentration (Cs-avg) at the outlet is given by: 


m 


Cone eer =. (18) 


where m is the total number of faces on the surface of the outlet and S, 


corresponds to the area with face index i, m7. 
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In the numerical simulation, the transient radon exhalation rate for diffusion was 


calculated as follows: 


E,_,=nDVC 


n 


x=0° (19) 


The transient radon exhalation rate for the convection is calculated by: 


E =nve 


x=0° (20) 


where v is the physical velocity of gas convection at the porous medium-air 


interface, ms". 
The total of the transient radon exhalation rate is defined by: 
E, = Epa tage (21) 


The initial radon exhalation rate, corresponding to the radon exhalation rate of 
the surface of porous medium in the natural state, is the value intent to be measured 
experimentally. In the numerical simulation of the homogeneous porous medium, the 
initial radon exhalation rate was equivalent to the value calculated using Eq. . 


The nonuniformity of the radon concentration distribution in the chamber is 


defined as?! 


NUI — Cras = Evin (22) 


v-avg 


where C-max is the maximum radon concentration in the chamber, Bq m”? and C,- 


min is the minimum radon concentration in the chamber, Bq m°. 


A lower NUI indicates better uniformity and vice versa. 


5 Results and discussion 
5.1 Effect of vent tube structure on radon concentration 

Prior to the experimental simulation, the vent tube structures of the chamber 
were designed and reconstructed. Two schemes were established, the vent tube with 
end opening for Case 1 and the vent tube with four side wall openings (82mm) on the 
wall for Case 2 (both under 1 L mint, 1x10° m°). The growth curves of the radon 


concentrations for these two cases are shown in Fig. 6a. In addition, the permeability 
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in Case 2 was changed to 1x10"! m? (Case 3 with four openings) to test the effect of 


permeability. 
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Fig. 6. (a) Growth curves of Cs-avg and C\-avg in the first 600 s (case 1 and case 2), (b) the full growth curve of radon 
concentration for the R model, and (c) exhibition of all sorts of radon exhalation rate. 


Prior to reconstruction (case 1), Cs-ave was larger than C)-avg in the first 100 s, 
with the gap narrowing and remaining for the rest of the time. However, Cy-avg was 
reckoned to grow faster than Cs-avg in view of the initial radon concentration being 
equal to 0 Bq m”, implying defects in the vent tube structure in Case 1. Upon 
modification (Case 2), Cs-avg increased more slowly than C,-avg and the discrepancy 
remained constant. C\-ag, representing the actual trend of radon concentration in the 


chamber, was used for subsequent simulations to study the migration of radon in the 


chamber. 
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Fig. 7. The streamlines and distribution of the radon concentration along the diagonal section of the accumulation 
chamber at 600 s for (a) Case 1, (b) Case 2, and (c) Case 3. 


The radon concentration contours for Cases 1 and 2 at 600 s are shown in Fig. 7. 
The streamlines in Case 1 show that the jet from the vent tube enters the porous 
medium and then returns to the chamber, directly affecting the radon concentration 
distribution in the porous medium, thus leading to the inhomogeneity of the radon 
concentration in the chamber. The streamlines in Case 2 mainly staying in the 
chamber indicate few impacts on the radon concentration distribution in the porous 
medium, in which the maximum, minimum, and average radon concentration in the 
chamber for Case 2 is 8676.49, 6208.08, and 7089.7 Bq m”, respectively, while those 
for Case 1 are 9473.19, 4444.73, and 6575.9 Bq m°. The NUI for Cases 2 and 1 were 
34% and 76%, respectively, revealing that the uniformity of the radon concentration 


distribution in case 2 significantly improved. 


5.2 Permeability effect on radon exhalation rate 

The £;,, on the surface covered by the chamber at 600, 1200, and 1800 s were 
1.03988, 0.84632, and 0.69968 Bq m?s! for Case 2, and 1.03943, 0.84595, and 
0.69937 Bq m?s! for Case 3. The approximate similarity of the radon exhalation rate 
as well as the contours and streamlines for Cases 2 and 3 imply that the effect of 
permeability on the closed-loop measurement of the radon exhalation rate can be 
neglected. Therefore, the permeability is set to 1x10! m? in the subsequent 


simulations. 
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5.3 Reference radon exhalation rate 

For assessment of radon exhalation rate under AC model, reference exhalation 
rate (Eye) under R model was calculated. The C,-avg growth curve from 0 to 1.2x107 s 
was obtained by radon accumulation until saturation in the sealed container (Fig. 6b). 
The initial (Ea), fitted (E¢e»), and total transient radon exhalation rates (En) were 
calculated, as shown in Fig. 6c. As the transient radon exhalation rate for convection 
remains at 0 Bq m? s! in approximation, the effect of convection on the exhalation 
rate measurement can be neglected. Thus, the total transient radon exhalation rate was 
determined by the transient radon exhalation rate for diffusion (En-a), which is 


represented by the overlap of En and En-a in Fig. 6c. 


The fitted radon exhalation rate characterized the average exhalation for a given 
duration. Upon saturation for radon accumulation, Ee, was approximately equal to 
En, suggesting the consistency of the fit, as well as the back-diffusion phenomenon, 
with Eyexp (0.520835 Bq m? s!) being only 35.59% equivalent to Ea (1.463235 Bq m? 
s!). 


5.4 Flow rate and insertion depth effect on radon exhalation rate 
The increased insertion depth of the chamber has been experimentally proven to 
reduce leakage during measurements. The effects of insertion depth (H1/H3/H5) and 
flow rate (0.5/1/2 L min`!) on the radon concentration and exhalation rate are analyzed 
in this section. The radon distribution in the accumulation chamber at 1800 s is shown 
in Fig. 8. At a low flow rate (0.5 L min‘), the radon concentration distribution is 
subject to notable spatial variation, while it gradually becomes uniform as the flow 
rate increases. In addition, lateral diffusion is responsible for the average 
concentration in the chamber for H1 being lower than that for H5, implying that the 


effect of insertion depth on the measurement cannot be neglected. 
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Fig. 8. Radon concentration distribution in the accumulation chamber at 1800s (flow rate with 0.5/1/2 L min“ and 
insertion depth with 1/3/5 cm). 


Excluding the points in the initial 60 s, the average value for every 60 s was 
calculated as the concentration at the midpoint, and further fitted to obtain the 
exhalation rate. For example, the average radon concentration from 61 s to 120 s was 
calculated as the value at 90 s. The duration for curve fitting included 600 s, 720 s, 
900 s, 1020 s, 1200 s, 1500 s, and 1800 s, while it was 240 s, 360 s, 480 s, 600 s, 720 
s, 840 s, 960 s, 1080 s, 1200 s, 1320 s, 1440 s, and 1560 s for linear fitting. 


Efexp and Epin were obtained by fitting the points to the curve and linearly fitting 
the equations using different time ranges, as shown in Fig. 9. A long duration (1800 s) 
is unfavorable for curve fitting, the discrepancy of which with E, enlarged (0.9% to 
3.3% increase in deviation for various insertion depths compared with 600 s). In 
contrast, the linear fitting value declined significantly with increasing duration (14.7% 


to 22.8% increase in deviation for various insertion depths from 240 s to 1800 s). 
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Fig. 9. The radon exhalation rate for the fitting value (Efexp in curve fitting and Efn in linear fitting) under various 
duration, the transient value in numerical simulation (En), and the analytical solution (E4) (flow rate with 0.5/1/2 L 
min", insertion depth with 1/3/5 cm, and permeability with 1x107!° m°). 

Both Efexp and Efiin increased with the flow rate for a fixed insertion depth closer 
to Ea. Throughout the accumulation duration, En maintained a decreasing trend and 
was only 48—59% (for various insertion depths) of E, at 1800 s, indicating that Ey is an 
averaged value for a certain duration and is unable to represent En at a specific 


moment. Moreover, Ey was more susceptible to back-diffusion for longer fitting 


durations. 


For an in-depth analysis of the insertion depth effect (1 L min"), Eye and Evin, 
as well as the effective decay constant (Ae), were calculated, shown in Fig. 10. A. 
(2.23582x10° s1) is comparable to Arn during closed-loop measurement in R model, 
in view of the assumption of radon-impermeable wall of the sealed container. A 
significant decrease of 27.2%-—34.2% is shown for 4e with the increased insertion 
depth, implying that the leakage was effectively reduced in line with the previous 


experimental observations.’ 
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Fig. 10. (a) Effective decay constant (4e) for H1, H3, H5, and R model using curve fitting (1 L min”, Arn=2.1x10-€ 
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Furthermore, the ratios of Ef and Ey to Ea were calculated (Fig. 10b and 10c). 
The ratio for the curve fitting (E~exp) is found to be positively correlated with 4e. 
Increased insertion depth results in lower 4e and the ratios, suggesting an obvious 
back-diffusion effect. Unlike curve fitting, the ratios for linear fitting showed slight 
growth. Specifically, the exhalation rates of curve fitting using 1800 s and linear 
fitting using 480 s for H1 were 0.89952 and 0.81364 Bq m? s! respectively, while 
they were 0.89241 and 0.83551 Bq m?s for H5. The effect of increased insertion 
depth on lower leakage is demonstrated by linear fitting (1.e., the fitted exhalation rate 
is closer to E, obtained by the increased insertion depth), whereas curve fitting for 


longer durations is more affected by back-diffusion. 


In general, a uniform mixture of radon inside the chamber facilitates the radon 
exhalation rate measurement as the flow rate increased from 0.5 to 2 L min’. 
Although the increased insertion depth helps lower the leakage rate, the side effect is 


the decreased accuracy of the curve fitting owing to back diffusion. 


5.5 Validation of the numerical model 
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The reliability and accuracy of the numerical model are validated. In previous 
experiments, Gutiérrez-Alvarez’!! deployed two reference boxes (RB1 and RB2) to 
obtain reference radon exhalation rate. Cylindrical (CX, where X denotes the insertion 
depth in centimeters) and rectangular (referred to as VX, where X has the same 
definition) accumulation chambers were designed to evaluate the influence of 
insertion depth. V10, C0, C3, C6, C9 were tested on RB1, whereas V10, C0, C3, C6, 
C9, C12, and C14 were tested on RB2. 


The effective decay constant measured on RB1 (approximately 4.8810 s-!) is 
62% greater than that on RB2 (approximately 3x10° s!) measured by RAD7, whereas 
it is 58% greater by Alphaguard, implying that the container of RB2 is better sealed 
during measurement. The effective decay constant for the accumulation chamber 
peaked at an insertion depth of 0 cm (CO), approximately 1.88<10~ s! on RB1 and 
2.5x10* s! on RB2. As the insertion depth increased, the effective decay constant 
declined by 38%—52% to almost 2.9x 10° s", indicating a lower leakage rate. The 
variation in the effective decay constant for the reference boxes (RB1 and RB2) and 
the accumulation chamber fit well with the simulated results shown in Fig. 10a. The 
effective decay constant decreased by 27.2%-—34.2% with increasing depth from 1 to 5 


cm. 


Furthermore, the ratios of the curve-fitted radon exhalation rate to the reference 
values for CO, C3, and C6 on RB2 were chosen for comparison with those for H1, H3, 
and H5 because the airtightness of RB2 was better than that of RB1. The ratios for CO, 
C3, and C6 were 1.65, 1.22, 1.1 measured by Alphaguard, and 1.18, 1.04, 1.02 by 
RAD7, respectively. An evident downward trend was observed for both Alphaguard 
and RAD7, which agrees well with the trends of the ratios for H1, H3, and H5 by 
exponential fitting in Fig. 10b. 


Although the reference radon exhalation rate (equivalent to the value obtained 
from the R model in this study) was used as an intrinsic value for comparison, the 
potential discrepancy caused by back diffusion should be noted. Inaccurate reference 


values can potentially render comparisons unreliable. 


In summary, the simulation not only validates the variation in the experimental 


data but also further reveals the positive correlation between the back-diffusion 
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phenomenon and the increase in insertion depth. 


5.6 Diffusion coefficient effect on radon exhalation rate 

In view of the inevitable back-diffusion phenomenon mentioned in Section 5.4 
(i.e., the curve-fitted radon exhalation rate declines with an increase in the insertion 
depth and with a lower effective decay constant), the effect of the radon diffusion 
coefficient for the porous medium on Ea, Ey, and En was additionally analyzed. 
According to Eq. , the radon exhalation rate is positively correlated with the effective 
diffusion coefficient, explaining why the Ea (1.2081 Bq m? s`!) for D1 scenario (the 
lowest effective diffusion coefficient among the four scenarios) is apparently lower 
than that for the other three scenarios (the former equivalent to 81.3%—87.7% of the 
latter in Fig. 11). Lower radon exhalation needs longer time for the radon 
concentration to reach saturated, evidenced by lower value of A. for D1 (Fig. 12a). 
Increasing the effective diffusion coefficient causes Ey or En to deviate from Ea, 
suggesting that the effective diffusion coefficient is an important factor in the back- 


diffusion phenomenon. 
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Fig. 11. The radon exhalation rate for the fitting value (Eye in curve fitting and Ezin in linear fitting) under various 
duration, the transient value in numerical simulation (£,), and the analytical solution (Ea) (D1/D2/D3/D4 
corresponding to 2.77 10°7/7.77*10°7/2.77*10°°/7.77*10°% m? s!, 1 L min", 1x10! m2, and 1 cm insertion depth). 

The ratios of Erto E4 using various fit durations for the four effective diffusion 
coefficient scenarios are shown in Fig. 12b and 12c. The ratios decrease with the 


increasing effective diffusion coefficient, going from 0.98 to 0.84 for curve fitting (Es 
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exp, 1800 s) and from 0.96 to 0.67 for linear fitting (Efun, 480 s). The deviations 
between the curve (1800 s) and linear (480 s) fitting were 2.08%, 3.88%, 9.55%, and 
20.68% for D1, D2, D3, and D4, respectively, making the exhalation rate for a lower 
effective diffusion coefficient that can be calculated approximately using linear fitting 


in a time-saving manner to replace curve fitting. 


108 | 600s 1200s MT 1800s 
| (a) 
1944 
2x10% f Arn 
10° 
D1 D2 D3 D4 
F E 600s B 1200s I 1800s i | 480s  :60s ER 1440s 


(b) (c) 


0.9 
0.8 a 
0.7 
0.6 
0.5 
0.0 
D1 D2 D3 D4 


Fig. 12. (a) Effective decay constant (4e) obtained by applying curve fitting for D1, D2, D3, and D4 scenario (1 L 
min", Arn=2.1x10° s!). Ratios of Esto Ea for D1, D2, D3, and D4 scenario using (b) curve fitting (E/exp) and (c) 
linear fitting (Epin) (1 L min`’). 
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6 Conclusion 
A CFD-based technique was employed to model the three-dimensional geometry 


and simulate various scenarios to measure the radon exhalation rate beyond the 
traditional method. The transient radon exhalation rate for diffusion maintained a 
decreasing trend during the initial 1800 s in the measurement owing to the back- 
diffusion caused by the accumulated radon, whereas that for convection contributed 
little to the radon exhalation. The results revealed that radon migration patterns are 
subject to the vent tube structure, effective radon diffusion coefficient, insertion depth, 
and flow rate. Their effects on the radon exhalation rate were analyzed, and the 


discrepancies between the initial, fitted, and transient exhalation rates were quantified. 
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The conclusions are as follows: 


The improved vent tube structure (tube with side wall openings) in the 
accumulation chamber was favorable for a stable radon concentration distribution in 
the porous medium, avoiding the disturbance of the jet from the tube with the end 
opening. The increase in circulation flow rate (from 0.5 to 2 L min`’) further improves 
the uniformity of radon concentration in the chamber (NUI from 76% to 34%), 
beneficial for accurate measurement of radon concentration. In addition, the 
permeability of the porous medium had little effect on the closed-loop measurement 


of radon exhalation rate. 


The effect of the effective diffusion coefficient on the radon exhalation rate 
measurement was notable. The deviation in the fitted radon exhalation rate (E) from 
the initial value (Ea) was small (max. 2.14%) for a relatively low effective diffusion 
coefficient (D1), whereas it increases for D4 (max. 15.96%), making the measurement 
susceptible to back diffusion. A lower effective diffusion coefficient indicates a lower 
radon diffusion capability in the porous medium, which limits the ability of radon to 
accumulate inside the chamber at high concentrations and diffuse into the porous 
medium. Therefore, the closed-loop method is more suitable for measuring the radon 
exhalation rate on the surface of a compact porous medium with a low effective radon 


diffusion coefficient. 


Although the increased insertion depth (from 1 cm to 5 cm) contributed to a 
decrease in the effective decay constant (a decrease of 23.17% and 19.55% 
corresponding to 600 s and 1800 s curve fitting, 1 L min‘'), back diffusion could not 
be avoided. This is evident from the deviation of the exhalation rate using curve 
fitting (Efexp) from the initial value (E4), which increased by approximately 7% (1800 
s), whereas that in the case of linear fitting (Ezin) decreased by approximately 12% 
(480 s). The deeper the insertion depth, the lower the leakage rate and, to some extent, 
the higher the equilibrium radon concentration inside the chamber. This, in turn, 
suppressed radon exhalation from the medium surface. Thus, linear fitting allows the 
radon exhalation rate to be calculated at a lower time cost and reduces the effect of 


back-diffusion. 
In summary, the consistency of the numerically simulated results with previous 
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experimental data signals that modeled the closed-loop measurement system of the 


radon exhalation rate by CFD-based numerical simulation is verifiable and reliable. 


The results are expected to provide theoretical guidance for the optimization of 


closed-loop measurement methods (e.g., circulation flow rate, insertion depth, and 


applicability of porous media). 


Nomenclature and units 


Ara Radium activity concentration, Bq kg"! 

C Radon concentration in the air, Bq m? 

Ca Radon concentration in the atmosphere, Bq m° 

Cs-avyg | Area-weighted average radon concentration at the surface of flow-out tube 
outlet, Bq m? 

C\-avg | Volume-weighted average radon concentration inside the accumulation 
chamber , Bq m° 

Cx-max | The maximum radon concentration value inside the accumulation chamber, 
Bq m? 

C\y-min | The minimum radon concentration value inside the accumulation chamber, 
Bq m°? 

Co Initial radon concentration, Bq m° 

D Radon diffusion coefficient in porous medium, m? s-! 

De Effective radon diffusion coefficient in porous medium, m°? s`! 

Din Molecular diffusion coefficient of radon in air, 1.05x 10- m? s`! 

Dp Average grain size, mm 

Dı Turbulent diffusion coefficient, m? s`! 

Di, Effective diffusion coefficient of porous medium for simulation calculation, 

D2, m? s! 

D3 

and 

D4 

E Radon exhalation rate, Bq m? s“! 

Ea Initial radon exhalation rate (equivalent to the analytical solution in 
diffusion), Bq m? s! 

Ef Fitted radon exhalation rate, Bq m? s! 

Efexp | Radon exhalation rate for exponential fitting, Bq m? s! 

Efn | Radon exhalation rate for linear fitting, Bq m? s! 

En Transient radon exhalation rate (En-a + En-c), Bq m? s! 

End Transient radon exhalation rate for diffusion, Bq m? s! 

Ene Transient radon exhalation rate for convection, Bq m? s- 

Ere¢ __| Reference radon exhalation rate, Bq m? s! 

Ern Radon emanation coefficient of porous medium, dimensionless 
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FVM | Finite Volume Method, a discretization method for numerically solving 
partial differential equations. 


g Gravitational acceleration, m s? 

Hl, The insertion depth of the accumulation chamber into porous medium, cm 

H3, 

and 

H5 

k Turbulent kinetic energy, m? s? 

K Permeability, m? 

NUI | None Uniformity Index, dimensionless 

R The reference model in which the reference radon exhalation rate is 
calculated. 


RB1, | The scenarios for measuring the reference radon exhalation rate. 35 kg 
RB2 | radium-containing medium with a height of 7 cm are deposited on RB1, 
while it is 70 kg with 14 cm on RB2. 


Sa The area of the surface where radon exhales, m? 


SRn Source term of governing equation of radon migration in numerical 
simulation, Bq m>? s! 


S: Area of the surface of the flow-out tube outlet, m? 

Scr Schmidt number, dimensionless 

t Time, s 

V Volume of the accumulation chamber, m? 

v Physical velocity of gas flowing in porous medium, m s`! 
x Vertical depth of porous medium, m 

a Free radon production rate, Bq m° s! 

E Turbulent dissipation rate, m? s? 

n Porosity, % 

ARn Radon decay constant, 2.1x10° s- 

Ab Back-diffusion rate, s`! 

he Effective decay constant, equivalent to ArntAvtu, s`! 

A The leakage rate for radon in the accumulation chamber, s-! 
u Dynamic viscosity, Pa s 

Lt Turbulent viscosity, Pa s 

ps Porous medium density, kg m? 

pa Air density, kg m° 

T Tortuosity factor, dimensionless 
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